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7.1 Equations of Motion in Rotating Coordinates 

Since the earth is rotating about the polar axis, the coordinate system fixed on earth is 
rotating. We need to know how to express the time rate of change of dynamical quantities 
in the rotating coordinates. 

A vector fixed in the rotating coordinate system is rotating in the fixed (inertial) coordi- 
nate system. Consider therefore a vector rotating in the inertial frame of reference. 

7.1.1 Vector of constant magnitude 




— # — * 

Figure 7.1.1: Vector A{t) rotating at the angular velocity Q. 

— * 

li A — AiCi has a constant magnitude but is rotating about an axis at the angular vecloty 
Q, what is the rate of change dA/dt in the fixed coordinate (inertial) system? Let 

dA = A{t + dt) - A{t) 

From Figure 7.1.1, 

^e\A\ sm 7— , 
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where subscript / signifies " inertial system" and e is the unit- vector along dA. Note e ± A 
and e _L Q. Hence, 

nxA 



and , 



Since 



dA 
Itt 



I Q X A 
VlxA 



VtxA 



dt ' 



^1 sm7^, 



Vtx A\=Vt \ A \ sin 7. 



it follows that 



dA 
Itt 



^VtxA. 



(7.1.1) 



In particular, let ^4 = Cj, i = 1, 2, 3 be a base vector in the rotating frame of reference, 

A^Ci^e 

Then 



dej 
dt 



Vt X e,-. 



(7.1.2) 



7.1.2 A vector of variable magnitude 

Let 

be any non-constant vector in the rotating frame, and let 

dB, ^ 



dB 
Itt 



R 



dt 



denote its rate of change in the rotating frame, then 



dB \ dB, 



dcj I dB 



dB 



dt 



dt 



ei + Bi—^=\ — \ + BiQx ei = \ — \ +nx B. 



dt 



dt 



R 



dt 



(7.1.3) 
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In particular, if S = r is the position of a fluid particle 

+ Qxf, 



df 
'di 



df 
'di 



R 
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Note that f is the same in any coordinate system. Now {df/dt)i is the velocity seen in the 
inertial frame of reference and {df/dt)R is the velocity seen in the rotating frame of reference, 
i.e., 

qi = qR + VLxr; (7.1.4) 

Next wc let q^i be the velocity vector of fluid in the rotating frame of reference; its rates 
of change in the two frames of reference are related by 



IR 



dt 



dqR 



dt 



+ QxqR. 



(7.1.5) 



R 



Taking the time derivative of (7.1.4), and assuming that the angular acceleration of earth to 
be zero, 

dn 



dt 



= 



we get 



dt 



' dqR 
dt 



' dq 



R 



^ dt 
'dqR 
dt 



+ fixqR + 0,x 



df^ 
dt 



+ Q xf 



2Q X (jR X {Q X f) 



(7.1.6) 



R 



Coriolis acc. centripetal 




Figure 7.1.2: Coriohs force, position vector and angular velocity 

The second term on the right is the Coriolis force, being perpendicular to both q and fl. 
The last term represents the centripetal force 

d X (O X f) = 

See Figure 7.1.2 for the geometric relations. 

The centripetal force may be written in terms of a centripetal force potential 0c where 

<f>c^^{nx f) • {n X r) = l\n\\-l. (7.1.7) 
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so that 



-V0, = -^e-l = -|Qprl, (7.1.8) 
ar I 



7.1.3 Summary of governing equations in rotating frame of refer- 
ence: 

Continuity: 

V-g = (7.1.9) 

In the coordinate system rotating at the constant angular velocity, the momentum equation 
reads, after dropping subscripts R 

p^^ + 2nx ^ = -Vp + pV(03 + (/)e) + l^V^q (7.1.10) 

where ^ 

7.1.4 Dimensionless parameters 



2flxq nU \flTj 

q- Vq Uyi U 



2nxq 2flL 



Rossby number 
Ekman number 



2VLxq 2W 2QL2 

= g 

For numerical estimate, we take Q = = 2.31 x lO^^s"^ and r = earth radius = 6400 
km. Then J^r ~ (2.31 x 10"^)^ x 6.4 x 10^ ~ 3 x lO'^m/s^ while g ~ lOm/s^ Hence 
g ^ Q^r; gravity is more important than centripetal force. 

7.1.5 Coriolis force 

Refering to the right of Figure 7.1.3 

Q = i'(-0 cos Q) + j(0) + kiSl sin Q) 

Introducing the spherical polar coordinates as in the left of Figure Refering to the right 
of Figure 7.1.3, with Q being the latitude. The Coriolis force is 
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Figure 7.1.3: The Northern hemisphere. 



2Q X g = 



-2Qcos^ 
u 



J 

o 

V 



k 

2n sin e 

w 



= I {-2flv sin 9)+ J {2nu sin 9 + 20w cos 9) + k (-20v cos 9) 

Consider shallow waters where the depth D is much less than the horizontal length L, 
i.e., D ^ L, and compare the two terms in the y direction of (j) 

2nwcos9 w ^ ^fD\ 
——— = -cot^ = O - cot^ < 1 
zSmsm^ u \LJ 

except along the equator where ^ = 
In the z direction of (A;), 

—2Q:VCOs9 —2Qucos9 D —2Qucos9 



IddP 
p dz 



Hence in shallow seas 
Define 

to be the Coriolis parameter, then 



2nxq'^ i{-2nvsm9) + j{2nusm9) 
f = 2Vt sin 9 



20, X —fv i + fuj 
In the northern hemisphere, < 9 < 7r/2. 



(7.1.11) 
(7.1.12) 



